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Abstract 

We consider a three spin component Bose Einstein Condensate as described 
by as many coupled nonlinear Schroedinger equations. For a very special 
ratio of the coupling constants, exact soliton solutions to this set of equa- 
tions are known. Here we find a simple representation including the N = 1 
solution based on the symmetry of the equations. This symmetry is described 
by a linear operator, the nonlinearity of NLS notwithstanding. Our useful 
representation opens the door to the nonintegrable case of general coupling 
constants. A new class of solutions is found. 



Bose Einstein condensates (BEC) offer an unusual opportunity to test 
quantum effects on the mesoscopic scale. In particular spinor BEC conden- 
sates furnish an excellent example of this possibility [1, 2, 3]. In 1998, the 
MIT group published observations of the first example of a spinor conden- 
sate, in which they transferred a spin-polarized ^^Na condensate, created in 
a traditional magnetic trap, into a dipole trap formed by the focus of a far- 
off-resonant laser [4]. Soon after that they observed the formation of spin 
domains resulting from population exchanges within the F = 1 manifold 
via spin-relaxation collisions. These collisions allowed for a redistribution 
of the spin population according to the constraints of the local magnetic 
field [5, 6, 7, 8]. Multiple condensates have numerous possible applications, 
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such as high resolution interferometry [9], multi channel signals and their 
switching, and even optodynamical devices [10] 

Here consider a spinor condensate formed by atoms with spin F = 1, 
described by a three components microscopic wave function. Some very 
special multi-component vector solitons have already been predicted [11, 12, 
13]. Their stability has been considered extensively, see for example [14, 
15, 16, 17]. In our previous work we investigated polar and ferromagnetic 
solitons, introduced in [11, 12, 13]. We showed that they can be generalized 
for a wider range of coupling parameters. By studying their collisions we also 
discovered a new class of solitary waves, which we called oscillatons. These 
include the polar and ferromagnetic solitons [18, 19]. In further studies we 
recognized that all bright solitons and solitary waves in spin F = 1 systems 
are oscillatons! This is discussed in the present work. Even though the 
collisions of oscillatons are inelastic, the emerging species are once again 
oscillatos [18, 19]. They form an "invariant subdivision". 

We consider a one dimensional, dilute gas of bosonic atoms with hy- 
perfine spin F = 1 at very low temperatures. Interactions between atoms 
can be parametrized by two coupling constants Cq = 47r^^(ao + 202) /3M 
and C2 = 4:7ih'^{ao — a2)/3M, where the a/ are s-wave scattering lengths of 
total spin / channels. The three component, macroscopic wave function, 
with m = —1,0,1, satisfies the Gross - Pitaevskii equation, which 
in dimensionless units has the form 

^^t^n^ = -ldl + J2\^rn'n<^m (l) 
\ m' / 

a \m',m" J m' 

Here /° (a = 1,2,3) are the angular momentum operators in a 3x3 repre- 
sentation and 7 = — C2/IC0I ^ 

This equation is symmetric under Galilean boosts and translations, t) = 



^If we choose as our length scale x n\co\ ^^'^ time scale t — > 2M ^^^j (7V 
is the number of particles and L is some characteristic length for given system, e.g. the 
transverse size of the trap confining a quasi one - dimensional condensate), all coefficients 
but the ratio between self and cross nonlinear coupling, — C2/|co|, which we denote by 7, 
will be equal to ±1. 
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B(j),Sx)^ = e'^^^^ — Sx — 2pt,t) and also under rotations, = 

e~^'^U{(j),9,ip)*^ = e~*'^e~*'^'^^e~*^-^^e~'"^-'^^^. Here, the operation U represents 
rotations of the quantization axis. In physical terms, U corresponds to turn- 
ing a Stern - Gerlach type apparatus used to measure the wave function 
components. We will put this symmetry to work. 

leda et al considered the case Cq = C2 = c < 0, equivalent to 7 = 1 [11, 12, 
13]. Equation (1) then is completely integrable. Bright A^-soliton solutions 
where found. For this case, binary collisions of solitons are illustrated and 
extremely long formulas describing them are given in [11, 12, 13]. The general 
one-soliton solution is given by 

2^e-*^* 

^sol 



1 + 4^e~^v^*^^"^"^ + Y^|(T)|^e^v^*^^~^°) 

(2) 



2y^e-v^(^-^o) X + ^— (f )*ev^(^-^'") fx 

where Xq = — ln(8/i/|(T)|)/2, x is a normalized polarization spinor x = 
(Xi,Xo,X-i) = EmlXml^ = 1) and f is a time reversing operator, 

defined as Tx = (X-d ^Xo? X+i)"^- The quantity (T) is a mean value of 
this operator: x^^X- O^^^ fi^^st goal was to find an expression for one-soliton 
solution, which is more transparent and expressed in terms of recognizable 
physical quantities. We arrived at: 

sech [y/m {x + Ax)] e''^+* 
$ = e-^^W(0,^,^) I I • (3) 

-y/fO sech [y/JH {x — Ax)] e*^~* 



Aside from rotation, this expression appears to be a superposition of two op- 
positely polarized displaced solitons. It reduces to leda's one soliton formula 
Eq. (2) when /i+ = /i_ = /i and Ax given by: 




Ax = -— In . (4) 



For a given value of (T) a direct relation between the angles of rotation and 
the components of the polarization spinor used in the papers of leda et al 
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Figure 1: Four examples of solutions to Eq. (6) for oscillaton components 77+ (continuous 
line) and tj^ (dashed line), 7 = ±1/2. 



[11, 12] is 

TT 1 - 

r = arer(T), 

2 2 ^ ' 

(3 = ^JJl^x — iip, 
Xo = e-'-^/\(^\ sin 9 cosh f3, 

For |(T)| —7- 1 and |(T)| — )■ single soliton Eq. (2) reduces to polar 
and ferromagnetic solitons^ respectively. In these limits, even without our 



■^In case of (T) — > the displacement Aa; tends to infinity and the spinor components 



X+i 
X-i 



V-2 

1 

7=2 



-i{t- 



m\ 



m\ 



cos I - I e 



^ — sin^ 



sm 



— cos^ 
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improvements Eq. (3) is recovered. In previous papers we only studied those 
two cases. The whole spectrum of one soliton solutions with < | (T) | < 
1 seemed to be fundamentally different then the above limiting cases and 
solutions of the form (3). This apparent distinction is now abolished. 

When Ax is arbitrary and /i+ is not equal to yU_ our Eq. (3) takes us 
out of leda's one soliton solutions as given by Eq. (2). To recreate this 
more general form in leda's terms, one would have to go to their two soliton 
solutions [11, 12]. 

Most importantly the form of Eq. (3) allows for an easy generalization to 
a much broader class of solutions when 7 7^ 1. This is a central result of our 
Letter. Its general form is 

y — ?7_(a;)e*^"* J 

For 7 7^ 1 the ?7's are in general no longer simple sech functions. Substi- 
tuting (5) into Eq. (1) yields the following system of two coupled ordinary 
differential equations, which the ?7's must satisfy (the prime denotes x differ- 
entiation): 

+ (1 + 7)r/| + (1 - l)vl + ^ = 0- (6) 

For each value of 7 there is a limit on possible values of the ratio of //_///+. 
The range of this ratio was found in [20] 

ro < ^ < - , for 7 > 

i < /£+ < ro , for 7 < ^ ' 

where 

In the limit 7 = 1, Eqs. (6) separate and the solutions are sech functions with 
arbitrary shifts and amplitudes, as mentioned above. In Fig. (1) we present 
four examples of rj pairs, for two different values of 7 = ±1/2 and zero shifts. 

In a numerical study we took a shifted solution given by (5) and slowly 
changed 7 starting with one. Even for small changes of 7 the components 



become independent. Each is a ferromagnetic soliton. 
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Figure 2: Schematic diagram of cases considered in this Letter. 

?7+ and rj^ attract or repel each other, due to the nonhnear couphng, and 
they begin to move with respect to each other. On top of that, we also 
observed inevitable excitations and radiation. However, if there are no shifts, 
an adiabatic change of 7 leads to adiabatic deformation of the ?7's and /I's, 
preserving the form of Eq. (5). 

Equation (6) is considered in the literature in the context of solitary waves 
in a two component system [20]. However, this is a merely formal similarity. 
In [20] the author considers a two component system, not supporting the 
possibility of a transfer between components. Here we focus on the F = 1 
manifold of a spinor condensate and this is a three component system. In 
our case the operator h( leads not only to the usual shift in population, but 
also to population oscillations between components. They are possible due 
to the existence of source terms in the Eq. (1), proportional to 7. Due to 
these oscillations we call these entities oscillatons (The name is not entirely 
new - oscillatons appear as solutions to the Einsten-Klein-Gordon equations, 
see E. Seidel et al, Phys. Rev. Lett. 66, 1659 (1991) and 72, 2516, (1994)). 
The explicit time dependence of the density of the three components of the 
oscillatons are given by 

|$o|^ = ^sin^e (r/^ + r/! + 2r/+r/„cos[a;t - 2?/']) , (9) 

9 9 1 

\^±i\'^ = sin^ - + cos'' - — -?7+?7_ sin^ ^ cos [w t — 2'?/'] , (10) 

which exhibit oscillatory behavior with frequency u = fi^ — fi^. Fig. 2 places 
the leda at al solution in a wider context of our family of calculations. For 
more details of oscillaton collisions consult our previous papers [19, 18]. 

Another case of interest is when 7 = 0. In many experimental applications 
7 is indeed very small (it is proportional to the difference between the two 
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characteristic scattering lengths Oq and 02). Take for example ^^Rb with 
7 ~ 0.01 [3]. When 7 is exactly zero, our system is once again integrable. 
There is no transfer between components of the wavefunction and the rotation 
U becomes a trivial operation. Notice that now tq must be exactly one 
and hence /i+ = There is no room for oscillations since the frequency 
= /i+ — yU_ is zero. This is not surprising, since for 7 = the equation of 
motion no longer admits spin mixing. 

In conclusion, when treating spinor Gross-Pitaevskii equations {F = 1), 
it has been possible to improve on the leda et. al. formula for a one soliton 
solution. The adopted form, apart from being simpler and easier to manipu- 
late, covers more cases and leads us into the world of oscillatons, entities with 
interesting properties including population oscillations between spin compo- 
nents following from the rotational symmetry of the system. These entities 
have been found recently [18, 19]. Now they follow from a systematic proce- 
dure. 
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